Fission properties for r-process nuclei 
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We present a systematics of fission barriers and fission lifetimes for the whole landscape of super- 
heavy elements (SHE), i.e. nuclei with Z > 100. The fission lifetimes are also compared with the 
a-decay half-lives. The survey is based on a self-consistent description in terms of the Skyrme- 
Hartree-Fock (SHF) approach. Results for various different SHF parameterizations are compared to 
explore the robustness of the predictions. The fission path is computed by quadrupole constrained 
SHF. The computation of fission lifetimes takes care of the crucial ingredients of the large-amplitude 
collective dynamics along the fission path, as self-consistent collective mass and proper quantum 
corrections. We discuss the different topologies of fission landscapes which occur in the realm of 
SHE (symmetric versus asymmetric fission, regions of triaxial fission, bi-modal fission, and the 
impact of asymmetric ground states) . The explored region is extended deep into the regime of very 
neutron-rich isotopes as they are expected to be produced in the astrophysical r process. 

PACS numbers: 21.60.Jz, 23.60.+e, 25.85. Ca, 26.30.Hj, 26.50.+X, 27.90.+b 
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I. INTRODUCTION 



The existence of super-heavy elements (SHE) above 
the naturally existing ones has attracted much attention 
in the past decades [l|, @] . The topic remains of high ac- 
tual interest as the new and heavier synthesized elements 
are added every year to the list, for a few examples from 
the rich list see 0-0] ■ Super-heavy elements are also pro- 
duced during the r-process [Sj-flQI] and their properties 
are important in order to determine the upper end of 
the nucleosynthesis flow. The key question in the study 
of SHE is their stability against the various decay chan- 
nels as a-decay, /3-decay, and particularly spontaneous 
fission. This paper aims at a theoretical survey of fission 
lifetimes for SHE. It will establish a systematics all over 
the landscape of SHE from the experimentally accessible 
neutron poor ones to the very neutron rich species which 
may occur in the r process. Fission lifetimes will also 
be compared with the lifetimes for a-decay. This survey 
is based on a theoretical description at level of a self- 
consistent mean field (SCMF). Such models came into 
practice about 40 years ago and have been steadily devel- 
oped to deliver now a reliable description of nuclear struc- 
ture and dynamics, for recent reviews see (ToT - H^ . We 
use here in particular the Skyrme-Hartree-Fock (SHF) 
approach which stays in the non-relativistic domain and 
employs an effective en ergy functional corresponding to 
zero-range interactions [12| . 

The first theoretical estimates of fission stability were 
performed on the grounds of the shell correction energy 
within phenomenologically adjusted shell model poten- 
tials [15( and studies of shell structure persist to be of 
high exploratory value also for self-consistent approaches 



[HI, [T3] ■ In the realm of SHE, one finds broad islands of 
shell stabilization rather than the narrow and deep val- 
leys as they are typically found for lighter nuclei fl8| . 
The emergence of large regions of stable nuclei is, in 
fact, favorable for the potential experimental accessibil- 
ity. The next step after estimates from the shell correc- 
tion energy is to check the systematics of fission barri- 
ers in SHE. There exists a wealth of information about 
fission in actinide nuclei (l9j which helps to probe the 
predictive value of the theoretical approaches, e.g. of 
the SHF method [2(| [H|. (In fact, a fission barrier was 
used in the calibration of one SHF functional [22j.) How- 



ever, comparison with experimental data requires to go 
beyond a pure SHF description and to take into account 
collective correlations (from rotation and low-energy vi- 
bration) which can modify the fission barriers by up to 
2 MeV [HI, Hij]. The systematics of fission barriers in 
SHE is simplified by the fact that there is only one fis- 
sion barrier to be considered (as opposed to actinides 
with their double- humped barrier) . It was found [2l], [25| 
that SHF provides estimates of islands of fission stability 
which are qualitatively in accordance with experiments 
[26U28I . The ultimate goal is, of course, to estimate the 
fission lifetimes directly. However, self-consistent calcu- 
lations of fission lifetimes are extremely demanding and 
thus have come up only recently, see e.g. f2p| - f3lT ] (which 
are mostly using still approximate masses and quantum 
corrections [23|) or [32j for a fully self-consistent calcula- 
tion. In this paper, we employ the method as presented 
in : 32| for establishing the systematics of fission lifetimes 
for all conceivable SHE. We will discuss the influence of 
the choice of the SHF parameterization and we will com- 
pute and compare also the lifetimes for a decay. 
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The paper is outlined as follows: Section |TT] presents 
the formal framework for the computation of fission life- 
times. Section [TTT] discusses a variety of results on fission 
barriers, fission lifetimes and a comparison with a-decay 
lifetimes. 



II. FORMAL FRAMEWORK 
A. SHF approach and pairing 

The basis of the description is SHF augmented by BCS 
pairing. The method is widely used and well documented 
in the literature, for reviews see e.g. [HI, Hj- We report 
here briefly the actual input and usage. 

The mean-field state is a BCS state characterized by a 
set of single-particle wavefunctions {(p a ,a = 1...0} and 
corresponding BCS occupation amplitudes v a . The SHF 
energy functional depends on density p, kinetic energy 
density r, spin-orbit density J, current j, spin density 
<r, and spin kinetic density t. We use the standard 
form of the Skyrmc functional, in some cases augmented 
by an isovector spin-orbit term [33|. We employ consis- 
tently the full form including the time odd currents (J, 
er, t) which play a crucial role in computing the col- 
lective masses along the fission path (see section IIIBI) . 
For the pairing functional we use the density dependent 
zero-range pairing force [34ll35j in the stabilized form (36j 
which reads in detail 
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[121 [T3|. Steadily growing availability of data from ex- 
otic nuclei and different preferences in the choice of the 
fit data has led to a variety of SHF parameterizations. 
One needs to check the results for sufficiently different 
parameterizations in order to explore the predictive value 
of SHF calculations. We will consider the following pa- 
rameterizations: SkM* [22j as a traditional benchmark 
because it is one of the first parameterizations deliver- 
ing a quality description of nuclear ground states; SkI3 
[33| which for the first time exploits the freedom of an 
isovector spin-orbit coupling thus simulating in this re- 
spect the situation in relativistic mean-field models; SLy6 
[40( which had been developed with a bias to neutron 
rich nuclei and neutron matter aiming at astrophysical 
applications; SV-min 41| as a recent development us- 
ing a large pool of semi-magic nuclei which were checked 
to have negligible correlation effects; and SV-bas which 
was adjusted to the same data as SV-min with an addi- 
tional constraint on nuclear matter properties (symmetry 
energy, isoscalar and isovector effective masses) to tune 
giant resonances together with ground state properties 
[41j . We will also show a result for HFB-14 as one rep- 
resentative in a large series of parameterizations derived 
in large scale fits biased on a comprehensive description 
of binding systematica, in this case referring to published 
data HI]. 

All these parameterizations employed different pairing 
recipes in their original definition. In order to make 
calculations better comparable, we use the same pair- 
ing functional ([1]) for all parameterizations and tune the 
pairing parameters (V P air,p, Vpair.n, Po, P air) to the data 



from even-odd staggering as summarized in [411 ] . Such a 
separate adjustment for each force is crucial because the 
actual pairing gaps energy from an interplay of pairing 
strengths Vp a ir,q and level density which depends sensi- 
tively on the effective mass and thus varies dramatically 
with the parameterization [l2|, [l4| . 



where u a = y/l — is the pair density and f a is a cut- 
off weight as defined in [3?1 ■ When proceeding along the 
deformation path, the pairing energy, impair, is plagued 
by a possible phase transition to the breakdown of pair- 
ing which, in turn, leads to singularities in the collective 
mass. A widely used method to stabilize pairing against 
breakdown is to employ the Lipkin-Nogami recipe, see 
e.g. [HI]. This, however, is not always effective enough 
in the breakdown regime. The stabilized functional (fTaj) 
provides a much more robust scheme. We use it with 
E cu tp = 0.6 MeV which leaves ground state properties 
nearly unchanged and is at the same time very efficient 
in suppressing singularities in the collective mass. 

The form of the SHF functional is more or less pre- 
scribed by a low-momentum expansion of a fictitious 
underlying effective two-body interaction [39| . But the 
model parameters cannot yet be derived from ab-initio 
methods with sufficient precision. It is customary then 
to adjust the parameters of the SHF functional to ex- 
perimental data, mostly from ground state properties 



B. Microscopic computation of fission lifetimes 

Fission represents a substantial rearrangement of a nu- 
cleus from one into two fragments. SCMF models are 
well suited to track this process in a least prejudiced 
manner. They require only one constraint to force a 
stretching of the system along the various stages while 
all other details of the rearrangements and shapes result 
automatically of the calculation. Actually, one uses an 
isoscalar quadrupolc constraint because the first stages of 
fission develop out of large-amplitude quadrupolc modes 
of the mother nucleus. Even this last piece of guess- 
work could be eliminated by using the recipes of adi- 
abatic time-dependent Hartree-Fock (ATDHF) [Hill]. 
This, however, has not yet been accomplished for the 
very heavy systems considered here. We stay at the level 
of constrained SHF and employ ATDHF only for a self- 
consistent evaluation of the collective mass and quantum 
corrections 12311 . 
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FIG. 1. (Color online) Upper panel: Raw potential energy 
surface (PES) V(a 2 o) and ZPE corrected PES V(a 2 o) for 
290 Sg. The energy Eo of the collective ground state is in- 
dicated by a heavy horizontal line. The points a, b, and 
c indicate the three crossing points of the Eq line with the 
collective potential V{a.2o). Tunneling and repetition rates 
W and T are also indicated. Lower: Inverse collective mass 
B a2Q calculated by a self-consistent cranking scheme (ATDHF 
cranking). Schematic nuclear shapes are sketched. 



Fission barriers have been discussed already in the 
early stages of SCMF models and have even been used as 
benchmark for calibration j22j . The calculation of fission 
lifetimes are much more involved as their computation 
requires not only the potential energy surface along the 
fission path, but also the corresponding collective masses 
and a safe estimate of the collective ground state corre- 
lations for the initial state. Thus, the vast majority of 
calculations of fission lifetimes employ the microscopic- 
macroscopic method which combines shell corrections 
with a macroscopic liquid-drop model background, see 
e.g. [HI, EH- Self-consistent calculations of fission life- 
times are still rare, see e.g. [29|, [3(|, and mostly use still 
approximate masses and quantum corrections. On the 
other hand, just because the computation of fission life- 
times is so demanding they serve as extremely critical 
observables probing all aspects of the effective nuclear 
interaction, its global bulk properties as well as details of 
shell structure. 

To calculate fission lifetimes within SHF we use the 
scheme as presented in (32[. We consider a fission 
path which evolves along axially symmetric shapes and 
characterize the deformation by the dimensionless axial 
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The steps to compute the fission lifetimes can then be 
summarized as: 

1. The fission path is a set of mean-field states 
{|<i> Q;20 )} representing the stages on the way from 
the ground state to fission. It is generated by 
quadrupole-constrained SHF (complementing the 
mean-field Hamiltonian by a constraining poten- 
tial, i.e., h — > h — A • Q2o)- 

2. The energy expectation value corresponding to 
|3> Q20 ) yields a "raw" collective energy surface, 
V(a 2 o)- 

3. The collective mass, B(a 2 o), an d moments of iner- 
tia are computed by self-consistent cranking (often 
called ATDHF cranking) along the states |$a 20 ) of 
the path [4f|. 

4. Approximate projection onto angular momentum 
zero is performed using the moments of inertia and 
angular-momentum width. 

5. Quantum corrections for the spurious vibra- 
tional zero-point energy (ZPE) are applied (using 
quadrupole mass and width). The result is the ZPE 
corrected potential energy surface (PES), V(a 2 o). 

6. The collective ground state energy Eq is computed 
fully quantum mechanically for the thus given col- 
lective Hamiltonian |46j |. 

7. The tunneling rate W at the given ground state 
energy and the repetition rates T are computed 
by the standard semi-classical formula (known as 
WKB approximation) using the quantum-corrected 
potential energy and collective mass (moments of 
inertia); the fission lifetime is finally composed from 
these two quantities as T/W: 



W = exp -2 
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T = H da 2 o ( \J B(a 2a )(E - V(a 20 )) 
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where the initial state would be classically bound 
in the interval (a, &), while the barrier extends over 
the interval (6, c) (see figure [1]). 

Point [5] in this list requires some explanation. The ax- 
ially symmetric fission path is described by three col- 
lective degrees of freedom (deformation a 2 o an d two ro- 
tation angles) while the full collective quadrupole dy- 
namics calls for the five-dimensional Bohr Hamiltonian. 
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FIG. 2. (Color online) Fission lifetimes of the isotopes No, 
Rf, Sg and Hs computed with three different Skyrme pa- 
rameterizations, as indicated, and compared with data from 
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In order to be consistent with the whole fission path, 
a three-dimensional collective dynamics for quadrupolc 
motion was derived which is restricted to axially sym- 
metric shapes. The method employs the norm and over- 
lap kernel of the topological Gaussian overlap approxi- 
mation [53, H| an d relies on the direct solution of the 
collective Schrodinger equation (see (1|), rather than es- 
tablishing a (reduced) Bohr Hamiltonian. 

Figure [TJ illustrates the collective parameter functions 
along the axially symmetric path for the super-heavy el- 
ement 290 Sg which are necessary for calculation of the 
fission half-lives. As can be seen in figure [T] the collective 
mass fluctuates strongly so it can hardly be approximated 
by some constant or weakly changing collective mass. 

As mentioned above, all calculations are performed in 
axial symmetry but allowing for reflection asymmetric 
shapes. This breaking of reflection symmetry becomes 
crucial in the outer region beyond the fission barrier. 
The ground states and the (first) barrier are usually as- 
sociated with reflection symmetric shapes with few ex- 
ceptions as discussed in the appendix. Breaking of axial 
symmetry towards triaxial shapes can occur in the barrier 
region. One knows from actinides that triaxial shapes can 
lower the barriers by about 0.5-2 MeV [Hill. Such low- 
ering is missing in axially symmetric calculations. The 
present results are thus to be understood as providing an 
upper limit for barriers and lifetimes. 



FIG. 3. (Color online) Fission barriers (lower) and life-times 
(upper) for two groups of experimental known super-heavy 
elements. Compared are results from a variety of Skyrme 
parameterizations with experimental data [H, l5oT - [58{ . The 
error bars on the barriers for SV-min are the uncertainties in 
the extrapolation as implied in the least-squares fits of the 
parameterization [4lj . 



III. RESULTS AND DISCUSSION 

A. Benchmark 

Lifetimes can be derived by the calculation scheme de- 
veloped above. Figure [2] shows fission lifetimes for four 
chains at the lower end of SHE. Three SHF parameteri- 
zations are compared. SkI3 which has a rather low effec- 
tive mass (m* I'm = 0.58) does not perform so well. This 
holds similarly for all forces with low effective masses. 
The results from the recent parameterizations SV-min 
and SV-bas which have effective mass m*/m = 0.95 and 
0.9 provide a more satisfying agreement taking into ac- 
count that an order of magnitude description is already a 
success for the extremely subtle observable of fission life- 
time. For the No isotopes, the agreement is acceptable 
in the average but there appears a strong deviation in 
the isotopic trend. This mismatch stems probably from 
the axial approximation and could well be explained by 
a strong isotopic change of the triaxial lowering of the 
barrier. 

Figure [3] shows experimental and calculated results on 
fission barriers (lower panel) and lifetimes (upper panel) 
for a few selected SHE, but now extending to heavier ele- 
ments and comparing more SHF parameterizations. The 
SHE represent two groups, one at the lower side (already 
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included in figure [2]) and another one with much heavier 
nuclei at the limits of present days data. The span of 
predictions from the various Skyrme forces is huge in all 
cases in spite of the fact that all these parameterizations 
provide a high-level description of nuclear ground state 
properties along the valley of stability. The variation 
of predictions may be a welcome feature as it provides 
additional selection criteria for a SHF parameterization. 
There remains, however, a problem when looking at the 
trend from the lighter side (Rf, Sg, Hs) to the heavier 
elements (Z=112, 114). All parameterizations produce 
a wrong trend of the predictions from the lower to the 
upper region. Barriers and lifetimes are well reproduced 
in the lower group by SV-min and SV-bas. But these 
parameterizations underestimate the barrier heights and 
lifetimes for the upper group [59| . The problem per- 
sists even with a more flexible density dependence of the 
Skyrme functional [o0| . It is also unlikely that triaxiality, 
ignored here, could help. It would worsen the situation 
for SV-min and SV-bas and the possible lowering about 
0.5-2 McV is insufficient to bridge the gap for the other 
parameterization. One has to keep in mind, however, 
that an experimental determination of lifetimes and bar- 
riers for the heaviest elements is a very demanding task 
and the data may not yet have reached their final stage 
such that the mismatch should presently not be over- 
interpreted. In any ease, we can expect from modern 
parameterizations as, e.g., SV-min a pertinent picture of 
the systematics of fission lifetimes for SHE. 



B. Fission topologies 



Panel c: shows a case where different fission paths 
emerge depending on whether one restricts the cal- 
culation on reflection symmetry or not or whether 
the calculation is going from oblate to prolate de- 
formation or vice versa. This suggests a struc- 
ture of two valleys where symmetric fission com- 
petes against asymmetric fission. This is also called 
multi- or bimodal fission and was already discussed 
in detail for self-consistent mean field models [Ifl 
[3ll . as well as for the microscopic-macroscopic 
finite-range liquid-drop model (63j (and citations 
therein). A full multi- modal treatment is presently 
beyond our possibilities. But the figure indicates 
that the barriers in the different channel are not so 
dramatically different. Thus it still provides a per- 
tinent picture if we consider one particular path, 
the one along asymmetric shapes starting from out- 
side. 

Panel d: shows the standard case which has one unique 
barrier in asymmetric calculations. In order to 
demonstrate the effect of asymmetry we compare 
with the PES from reflection symmetric calcula- 
tions. The latter show the double- humped fission 
barrier as it is known from actinides [64| . The al- 
lowance of asymmetric shapes removes the second 



barrier which holds for practically all SHE [21 



The overview demonstrates the large variety of topologies 
for the fission PES. This inhibits an automatic barrier 
search. Thus for the greater part of the investigated PES 
the minima and maxima relevant for the fission process 
were determined manually. 



When going through the variety of SHE, one en- 
counters much different shapes of the fission landscape. 
Four different situations can be distinguished as shown 
schematically in figure 0] 

Panel a: shows a case where the potential energy sur- 
face has a strong oblate minimum and optionally 
a secondary prolate minimum. It is thus assumed 
that the fission path is going through the triaxial 



plane [61 



In fact, we have seen from fully triaxial 
calculations that in most cases already the ground 
state acquires some triaxiality. Therefore, it is not 
possible to make reliable predictions using an axial 
code in such cases. No fission barriers and lifetimes 
will be shown for this region in the following sec- 
tions. 

Panel b: shows a PES for the case of a reflection asym- 
metric ground state which occurs in two region in 
88 < Z < 96 and 132 < N < 136 and again in 
96 < Z < 108 and 186 < N < 198. The breaking 
of reflection symmetry enhances the binding which, 
in turn, increases the fission barrier and thus leads 
to enhanced fission lifetime. This aspect will be 
discussed in more detail in the appendix. 



C. Systematics of fission barriers 

A simple first indicator of fission instability is the 
height of the fission barrier. Figure [5] shows the system- 
atics of fission barriers of SHE in the range 82 < Z < 120 
and 120 < N < 260. All elements are found to be stable 
against immediate nucleon emission at the ground state 
and along the whole fission path. In case of a double- 
humped structure of the barrier (commonly appearing 
for Z < 100) the higher barrier is plotted. Results are 
shown for the four different Skyrme parameterizations 
SLy6, SkI3, SV-min and SV-bas. 

At first glance it is apparent that the Skyrme forces 
SkI3 and SLy6 yield notoriously much higher barriers as 
SV-min and SV-bas. This can be traced back to a differ- 
ence in the effective mass m* /m. SLy6 and SkI3 have a 
very low mass m* /in = 0.69 and m* /in = 0.58, respec- 
tively while SV-bas and SV-min have effective masses 0.9 
and 0.95. A low effective mass leads to a too low density 
of single-particle states, and thus to larger shell correc- 
tion energies which, in turn, yield larger barriers. 

All forces provide the same trends over the landscape 
of SHE. There is a strong variation in fission barriers 
corresponding to the strong variations of shell structure 
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FIG. 4. (Color online) Schematic survey of fission modes for the various SHE. The two panels on the left and right show 
examples for the different regions as indicated: a) triaxial fission b) asymmetrical ground-state energy c) bi- or multi-modal 
fission d) asymmetric calculation removes completely the outer barrier. 
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FIG. 5. ( Color online) Systematics of fission barriers of even- 
even nuclei for four different Skyrme parameterizations SLy6 
[U, SkI3 [H, SV-min and SV-bas HO]. Grey indicates 
a possible triaxial fission, as already seen in figure [4] 



in the landscape of SHE. Several regions of high barri- 
ers occur. For low Z, i.e. Rn (Z=86), Ra (Z=88), and 
Th (Z=90), one sees two islands of high barriers, one at 
TV ~ 126 and another one at TV 184. Neutron numbers 



in between cover a region of lower barriers. They are 
particularly low for SkI3 and to some extend SLy6 while 
the fluctuations between high and low barriers are less 
dramatic for SV-bas and SV-min, corresponding to their 
generally smaller shell corrections. Stepping up to higher 
proton numbers Z there follow two more regions of high 
barriers, one of deformed SHE around Z/N = 104/152 
and one of spherical SHE around Z/N = 120/184 (spher- 
ical shell closure). The magic neutron number N=184 is 
clearly visible while an expected magic proton number 
near Z = 120 is indicated by a broad island of enhanced 
stability barriers around Z/N = 120/184 [H|. 

Considering the regime of nuclei relevant for r-process 
nucleosynthesis it is interesting to notice the appearance 
of a region of low fission barriers for Z ~ 84 in mov- 
ing from TV = 126 to TV = 184 for SLy6 and particu- 
larly for SkI3. These low fission barriers may allow for 
neutron-induced fission to occur as the nucleosynthesis 
flow moves from the N = 126 region to the TV = 184. 
The situation is different for SV-min and SV-mass, that 
predict much larger fission barriers in the r-process rele- 
vant region. For these parameterizations fission will only 
be relevant once the nucleosynthesis flow overcomes the 
TV = 184 magic number. All four forces agree in pre- 
dicting a rapid decrease of barrier heights going from the 
shell closure TV = 184 up to neutron-rich nuclei. This 
suggest a substantial decrease in the production of nuclei 
beyond TV = 184 during the r process. These aspects will 
be explored in a forthcoming publication. 

To complete the picture, figure [S] shows results from 
four other models for which data are publicly available. 
The calculation were performed by using the Extended 
Thomas- Fermi plus Strutinsky integral (ETFSI) [d5| . 
the Thomas-Fermi (TF) [66[ method , the macroscopic- 
microscopic finite range liquid-drop model (FRDM) [631 ] , 
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FIG. 7. (Color online) Systematics of fission lifetimes calcu- 
lated in the range of proton numbers 100 < Z < 120 and 
neutron numbers 140 < N < 260 for even-even nuclei us- 
ing the Skyrme parameterizations SkI3, SLy6, SV-min and 
SV-bas. 
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FIG. 6. (Color online) Results for fission barriers of four 
theoretical mass models HFB-14 g|, ETFSI [H, FRDM [||, 
TF in the same region of nuclei. The shown data are 
reduced to even-even nuclei to enable a better comparison. 



or the SHF approach with the parameterization HFB14 
[42l ]. All four theoretical mass models shows similar 
trends as observed in our studies. They confirm the 
region of high barriers below Uranium, around neutron 
number N = 184 and the island around Z/N = 104/152, 
the latter, however, less strongly developed in case of 
TF and FRDM. The three non self-consistent models 
(TF, ETFSI, FRDM) shift the third island of stabil- 
ity (Z/N = 120/184) down towards proton number 
N = 114. The rapid fall-off beyond N = 184 is con- 
firmed. The case HFB14 belongs to the SHF family. In 
spite of the narrow range of results one can conclude that 
its systematics is very similar to the other SHF cases. 
The actual barrier height are closer to SV-min and SV- 
bas (sometimes even below) . This is not surprising as all 
models in figure [5] have a large effective mass, equals or 
around m*/m = 1. 



D. Fission lifetimes 

Figure [7] shows fission lifetimes calculated with the 
recipe as outlined in section Hi Bl in the range 100 < Z < 
120 and 140 < N < 260 (from proton to neutron drip 
line) for the same Skyrme parameterizations as in the 
barrier systematics. Calculation of lifetimes were not per- 
formed for nuclei with proton numbers Z < 98, because 



the doubled humped barrier makes the evaluation of life- 
times cumbersome. The difference in barrier heights from 
to 12 McV (figure [5]) translates to a difference in life- 
times from almost immediate decay to 10 12 s and longer, 
demonstrating again the enormous sensitivity of fission 
lifetimes to all details of the model and computation. 

At first glance, the basic pattern resemble much the 
systematics of barriers. Long lived SHE are obviously 
found in the islands of high barriers. The island around 
Z/N = 104/152 is even broadened to higher Z and N 
towards the assumed neutron shell closure at N=162, 
especially for SV-min and SV-bas. This demonstrates 
that not only barrier but also barrier width and collec- 
tive mass can have a decisive influence. All parameteri- 
zations show a broad and deep valley of fission instability 
starting abrupt with neutron number TV = 186/188. If 
the r-process nucleosynthesis flow is able to overcome the 
N = 184 magic number, it will proceed by the region of 
large spontaneous fission lifetimes in figure [5j However, 
once the neutrons are exhausted and matter beta-decays, 
the region of short spontaneous fission lifetimes will be 
reached and no long live SHE will be produced. The sit- 
uation may be different for SkI3 and SLy6, depending on 
the extend of the region of short lifetimes above Z > 120. 

There is a large difference between the SHF parame- 
terizations in overall lifetime for elements with N < 190. 
The nuclei are much more stable for SkI3 and SLy6 than 
for SV-bas and SV-min. This is, of course, related to the 
overall difference in barrier heights (sec figure [5]) which 
can be traced back to different effective masses m*/m. 
This produces here even a qualitative difference: The 
parameterizations SV-min and SV-bas show a valley of 
fission instability between the islands around Z = 120, 
N = 180 and Z/N = 104/152 while SkI3 and SLy6 do 
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FIG. 8. (Color online) a-decay half-lives computed by using 
the Q a values and the semi empirical Viola Seaborg formula. 



not. The immediate consequence is that SkI3 and SLy6 
predict uninterrupted chains of a decay from the heaviest 
SHE down to actinides while SV-min and SV-bas have 
these a chains terminated by spontaneous fission. The 
latter is what is empirically found Q. The competition 
with a decay is discussed in section UlI El 



E. a-decay 

The a-decay half-lives are evaluated using the Viola- 
Seaborg relationship [67], IH| : 



log(r Q /s) = (aZ + b)(Q a /McY)- 1 ' 2 
+(cZ + d) + h\ os 



a = 1.66175, b 
c = -0.20228, d -- 



-0.5166, 
-33.9069, 



(5) 
(6) 

(7) 

(8) 



Q a = E(N—2, Z — 2) + £7(2, 2) — E(N, Z) , (9) 



E(2, 2) = £ cxp ( 4 He) = 28.3 MeV 



(10) 



and h\ og = for the even-even nuclei considered here. 
This requires as only input the Q a values which can be 
determined easily as difference of ground state binding 
energies. The latter are computed allowing for axial de- 
formations as well as reflection-symmetry breaking and 
including approximate angular momentum projection for 
deformed nuclei. It is to be noted that the difference of 
binding energies, as the Q a value, are predicted rather 
reliably with the SHF models although the binding en- 
ergies as such are notoriously underestimated for SHE 
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Figure [8] shows the systematics of a-decay half-lives 
calculated with the SHF parameterizations SLy6 and SV- 
min. In contrast to fission lifetimes, a-decay half-lives 
vary in general smoothly and steadily with a tendency 
to increase when going in direction of neutron rich SHE. 
An exception arc the spherical neutron shell closures at 
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FIG. 9. (Color online) Minimal lifetime plots for the compet- 
ing decay channels a-decay (green), and spontaneous fission 
(red). Experimental known decay channels [5£j are compared 
to results using the parameterization SV-min and SLy6. 



N=126 and N=184 which are clearly marked by a sudden 
decrease of a half-lives. But there is no detailed struc- 
ture like the islands of stable nuclei in spontaneous fis- 
sion systematics. Most of the nuclei in the shown region 
are very stable against a decay. It is only the band of 
neutron-deficient SHE at the left side of the region where 
a decay plays a role as competitor to fission and (3 decay. 
Comparing with the fission lifetimes in figure [7] that a 
decay prevails in any case for the islands of fission stabil- 
ity around Z = 120, TV = 180 and Z/N = 104/152. The 
parameterizations SV-min and SV-bas produce the pro- 
nounced valley of fission instability between these islands 
for which then fission takes the lead over a decay. This 
does not happen for SLy6 and SkI3 with their generally 
longer fission life times. 

A direct comparison of a decay and fission is provided 
in figure [9] showing the systematics of the dominant de- 
cay channels. Results are shown for the two parameteri- 
zations SV-min and SLy6 and compared with data. The 
experimental situation is just in between the two theoret- 
ical predictions. SLy6 produces too much fission stability 
thus giving a decay too much dominance while SV-min 
slightly underestimates the impact of a decay. Consider 
the a decay chain from Z = 118, TV = 176. It is ter- 
minated by fission already for Z = 116 for SV-min while 
the experimental chain terminates later at Z = 112. The 
results for SkI3 are very similar to those of SLy6 and the 
results for SV-bas to those of SV-min. 

It is to be noticed that the trends which are seen here in 
the systematics had already been spotted in figure[31 The 
parameterizations SV-min and SV-bas perform generally 
better, but are still plagued by producing a wrong trend 
when stepping to the heavier end of SHE where they yield 
too low barriers and lifetimes. A better compromise has 
yet to be worked out. 
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IV. CONCLUSIONS 

We have explored the systematics of fission barriers 
and lifetimes in the realm of super-heavy elements (SHE) 
on the grounds of self-consistent calculations using the 
Skyrme-Hartree-Fock (SHF) approach. The fission path 
has been generated with a quadrupole constraint pro- 
ducing a series of axially symmetric deformations while 
allowing for reflection-asymmetric shapes. The corre- 
sponding collective mass is computed by self-consistent 
cranking (often called ATDHF cranking). The quantum 
corrections to the collective potential (angular momen- 
tum projection, vibrational zero-point energy) are prop- 
erly taken into account. The fission life-time is com- 
puted for thus given potential and mass by the semi- 
classical WKB approximation, while the ground state en- 
ergy, which is at the same time the entrance energy for 
fission, is computed quantum mechanically in the given 
collective geometry of one axial deformation plus two ro- 
tation angles. Results have been produced for a couple 
of different SHF paramctcrizations to explore the sensi- 
tivity to the parameterization. For comparison, we have 
also computed the a-decay life-times using the Viola sys- 
tematics. 

A first test was performed by comparing with known 
fission life-times in isotopic chains in the lower region of 
super-heavy elements, Z = 104-108 and for the few avail- 
able data points in even heavier elements. The span of 
predictions is large whereby the effective mass of the un- 
derlying parameterization plays a crucial role. Satisfying 
agreement is found for modern parameterizations using 
effective mass around m*/m = 0.9-1. There remains, 
however, one open problem with the global trend: The 
parameterizations which perform almost perfectly in the 
region Z = 104-108 underestimate fission barriers and 
lifetimes in the heavier region Z = 112 and 114. 

The landscape of SHE separates into regions of differ- 
ent topology of the fission path. The most widely found 
standard case is a unique, axially symmetric fission path 
showing only one fission barrier; the second barrier which 
is known from actinides is suppressed by reflection asym- 
metric shapes which regularly develop for larger deforma- 
tions. Proton rich isotopes in 100 < Z < 108 show often 
a tendency to bi- or multi-modal fission where different 
fission and fusion paths compete. A small region around 
Z = 118 and N = 200 has oblate (if not triaxial) ground 
states and can decay only through a manifestly triaxial 
fission path. There are, furthermore, two small regions 
where the ground state is reflection asymmetric. This 
was shown to enhance barrier and lifetimes at a quanti- 
tative level, but not changing the global trends. 

The systematics of fission barriers and lifetimes shows 
the known islands of stability around Z/N = 104/152 
in the region of a deformed shell closure and around 
Z/N = 118/178 in a region of spherical isotopes. The 
actual values of barriers and lifetimes depend very much 
on the SHF parameterizations. Those with low effec- 
tive mass (here SkI3 and SLy6) produce very high bar- 



riers and lifetimes while those with high effective mass 
(SV-min and SV-bas) yield moderate barriers and life- 
times. The latter group also produces a valley of fission 
instability between the two islands, qualitatively in ac- 
cordance with the empirical findings. The way to very 
neutron rich ?*-proccss nuclei with N > 184 starts out 
with a large region of fission instability. Some stability 
is gained at the extremely neutron rich end. This makes 
unlikely the production of long-lived SHE above Z = 100 
by the r process. However, a more realistic estimate of 
the production of SHE by the r process will require nu- 
cleosynthesis calculations based on the present barriers 
and lifetimes. This will be the subject of a forthcoming 
publication. 

We have also computed a-decay lifetimes using the 
Viola-Seaborg formula. While the fission life-times show 
dramatic variation over the chart of super-heavy elements 
(from instability to Tfi ss = 10 16 s), the a-decay times 
vary gently with small overall changes and without vis- 
ible shell effects. The general crossover from a-decay to 
fission along the decay chains from the upper island of 
SHE is qualitatively reproduced by the family of SHF 
parameterizations with high effective mass. A quantita- 
tively reliable prediction of the switching point is detail 
on which the models have yet to be refined. 



Appendix: Symmetric vs. Asymmetric (incl. 
isomer) 

Reflection asymmetric shapes are a key issue in fission 
of SHE and thus have been much debated under different 
aspects as, e.g., suppression of the second barrier, im- 
pact of ground state asymmetry on the first barrier, or 
influence on bi- and multi-modal fission [3l|, [62|, [6^, [6t| . 
We will discuss here the effect of reflection asymmetric 
ground states on the fission barrier and subsequently on 
lifetimes. 

Figure \TU\ shows the PES of the two isotopes 286 Rf 
and 294 Rf. In the majority of the cases symmetric and 
asymmetric calculations provide the same ground state 
energy. This is illustrated by 286 Rf. In contrast, for 294 Rf 
the asymmetric calculation yields an energetically more 
favorable ground state (around a2o = 0.08), while in the 
symmetric PES it is difficult to locate the minimum (pos- 
sibly around ct2o = 0.32). It is to be remarked that the 
tendency to symmetry breaking is confined to the ground 
state region thus lowering the ground state energy. The 
absolute height of the barrier remains almost unaffected. 
As a consequence, an asymmetric ground state will lead 
to higher (relative) fission barriers. 

Besides Rf, the elements Pu, U, Th, R and Rn are 
also known for the importance of the octupole degree of 
freedom ([7(j and [lH). Figure Qj] summarizes for all 
relevant elements the differences of binding energies and 
fission barriers between reflection symmetric and asym- 
metric calculations for SLy6 (upper panels) and SV-min 
(lower panels). All these isotopes display basically the 
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FIG. 10. (Color online) PES of the isotopes 294 Rf and 294 Rf 
for symmetric and asymmetric shapes calculated with the 
Skyrme parameterization SLy6. 



effect as discussed for Rf, namely that the asymmetric 
shape affects predominantly the ground state thus lead- 
ing to an increase in the fission barrier which corresponds 
directly to the lowering of the ground state. The size 
of the effect changes quickly with proton and neutron 
number which is no surprise because symmetry break- 
ing is driven by (quickly changing) shell structure. This 
also explains that the results from SLy6 and SV-min are 
quantitatively so much different. SLy6 has a significantly 
lower effective mass than SV-min thus lower level den- 
sity and, in turn, larger shell corrections. The lowering 
of the ground-state energy was also investigated using 
other Skyrme forces [72j , where the results show a strong 
dependence on the effective mass m*/m and the pair- 
ing strength. A small effective mass or a small pairing 
strength lead to a big effect on the ground state energy 
and vice versa. 

Figure [TJ] compares the systematics of barriers (upper 
panels) and lifetimes (lower panels) with and without al- 
lowing for asymmetric shapes in the region of relevant 
isotopes. One spots a slightly increased fission stability 
in case of allowed asymmetry in the feature that the sta- 
ble regions are somewhat extended. However, the overall 
trends and the general impression of the plots of sys- 
tematics remains the same. The effect of asymmetric 
shapes appears rather at a quantitative level. The exam- 
ple demonstrates how robust the analysis of global trends 
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FIG. 11. (Color online) Difference of ground state energies 
AEo and fission barriers SBf between reflection symmetric 
and asymmetric calculations for the isotopes as indicated. 
The upper panels show results computed with the Skyrme 
parameterization SLy6 and the lower panels with SV-min. 
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FIG. 12. (Color online) Comparison of symmetric and asym- 
metric fission barriers and half-lives in region b) of the 
schematic survey figure [4] where asymmetric ground states 
can play a role. The upper panels show barriers and the lower 
panels lifetimes. The left panels show results from reflection 
symmetric calculations and the right panels from calculations 
where asymmetric shapes were allowed. The results are cal- 
culated with the parameterization SV-min. 
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